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Abstract: We perform a Next-to-Leading order analysis of AS = 2 processes beyond 
the Standard Model. Combining the recently computed NLO anomalous dimensions 
and the B parameters of the most general AS* = 2 effective Hamiltonian, we give 
an analytic formula for AMk and Ek in terms of the Wilson coefficients at the high 
energy scale. This expression can be used for any extension of the Standard Model 
with new heavy particles. Using this result, we consider gluino-mediated contributions 
to AS = 2 transitions in general SUSY models and provide an improved analysis of the 
constraints on off-diagonal mass terms between the hrst two generations of down-type 
squarks. Finally, we improve the constraints on R-violating couplings from AMk and 
Ek- 

Keywords: Supersymmetric Standard Model, Kaon Physics, CP violation, NLO 
Computations. 



1. Introduction 


The prescription of minimality in the number of new particles introduced to supersym- 
metrise the Standard Model (SM), together with the demand of conservation of baryon 
and lepton numbers, does not prevent the appearance of more than 100 new SUSY 
parameters in addition to the 18 already present in the SM. Fortunately for SUSY pre- 
dictivity, most of this enormous parameter space is already phenomenologically ruled 
out. Flavour changing neutral current (FCNC) and CP violating phenomena are the 
protagonists of this drastic reduction of the SUSY degrees of freedom. Therefore, the 
interest on the FCNC constraints on SUSY is well motivated: given the large degree of 
arbitrariness in the construction of low-energy SUSY models, they represent one of the 
main criteria we have for selecting viable theories. At the same time, they can shed 
some, albeit dim, light on the more fundamental physics from which low-energy SUSY 
stems. 

The ignorance of the underlying physics compels us to hnd ways to analyse the im¬ 
pact of FCNC and CP violating processes on low-energy SUSY in a model-independent 
manner. Since the pioneering work of Hall, Kostelecky and Raby in 1986 [1], the most 
adequate tool for such analyses has proven to be the so-called mass insertion approx¬ 
imation. A brief description of this method will be provided later on. We note here 
that this method has been applied, so far, mainly to the most genuine SUSY sources of 
FCNC, i.e. the Flavour Changing (FC) couplings of gluinos and neutralinos to fermions 
and sfermions [2]. Differently from the other sources of FCNC in SUSY (namely W, 
charged Higgs and chargino exchanges), these couplings have no analogue in the SM. 

As long as we deal with general squark mass matrices, the inclusion of gluino- 
mediated FCNC diagrams alone is sufficient to get the correct bulk of the SUSY con¬ 
tribution. In specihc models, corresponding to particularly restricted squark masses, it 
may occur, however, that other contributions become more important. This is the case, 
for example, of the Constrained Minimal Supersymmetric Standard Model (CMSSM) 
with exact universality of the soft breaking terms: in this case, in regions of the SUSY 
parameter space where charginos, stops and/or charged Higgs are relatively light, their 
exchange may overwhelm the gluino contribution to FCNC (see for instance [3]). In this 
work we intend to study the general case of squark mass matrices and, thus, we con¬ 
centrate on the gluino and squark exchange contributions. We postpone the inclusion 
of the chargino-squark contributions in the mass-insertion framework to a subsequent 
work. 

In the literature there exist four analyses which make use of the abovementioned 
method to perform a complete study of the whole variety of FCNC processes in the 
hadronic and leptonic sectors [3]-[6]. These studies do not include QCD radiative 
corrections and make use of the Vacuum Insertion Approximation (VIA) in the eval¬ 
uation of the hadronic matrix elements. A major step forward was recently made by 
the inclusion of the leading QCD corrections in the evaluation of the gluino-exchange 
contributions to mixing [7]. What Bagger et al. found is that the leading 



QCD corrections affect the results in a non-negligible way: in particular for compa¬ 
rable values of the squark and gluino masses, the QCD corrections increase the lower 
bound on the squark masses of the hrst two generations by a factor three. As for the 
S-parameters of the hadronic matrix elements, the VIA values were used. 

Motivated by the encouraging results of ref. [7], by the recent computation of the 
Next-to-Leading (NLO) QCD corrections to the most general AF = 2 effective Hamil¬ 
tonian [8] and by the lattice calculation of the full set of H-parameters con¬ 

tributing to the mixing matrix elements [9, 10], we present here the results of a 
new analysis of mixing which includes these progresses intended to improve the 

theoretical accuracy. 

A full NLO computation needs the evaluation of the 0{as) corrections to the Wilson 
coefficients at the scale of the SUSY masses running in the loop. This piece of infor¬ 
mation is not available yet for gluino contributions^. We can only invoke the smallness 
of as evaluated at such large scale as a (arguable) hint for the relative smallness of 
these corrections. Even with this limitation, we think that it is relevant to provide 
the improvement that the knowledge of the NLO QCD corrections to and the 

lattice values of the H-parameters allow us to obtain. 

A hnal comment is in order. One may wonder whether such a rehnement (inclusion 
of NLO QCD corrections and H-parameters from lattice computations) is worth in a 
study where many uncertainties, due to the presence of undetermined SUSY parameters 
and to the use of certain approximations (in particular the omission of the full set 
of FCNC contributions with possible non-negligible interference effects), are present. 
We think that the answer is positive at least for two reasons. The practical one is 
that, due to the complete model-independence of our approach, it is possible to “test” 
rapidly the impact of FCNC processes on new classes of low-energy SUSY models 
by just comparing their predictions for the FC parameters S (see below) with the 
constraints that we provide. Obviously a determination as precise as possible of these 
Ss is welcome. The second motivation is that with the future inclusion of the chargino 
exchange contributions to FCNC processes we will be able to obtain a comprehensive 
framework where to perform general studies of low-energy SUSY. In view of this, the 
use of consistently renormalised at the NLO and a better determination of the 

H-parameters play a major role in obtaining an efficient tool of analysis. Our aim is to 
approach as much as possible the situation of the computation of AF = 2 processes in 
the SM, where the complete LO [12] and NLO [13] analyses have been carried out. 

To this aim, we provide an analytic formula for the NLO in terms of the 

Wilson coefficients at the SUSY scale. This expression, together with the values of 
the H-parameters, can be readily used to compute AMk and Ek in any given SUSY 
model. Indeed, the formula we give is valid in any extension of the SM in which new 
heavy particles contributing to FCNC processes are present, since the basis of operators 

^NLO QCD corrections to the matching conditions for charged-Higgs and chargino contributions 
to £k have been very recently computed [11]. 



considered here is the most general for AS = 2 transitions. One just has to plug in 
the initial conditions computed in his favourite model. At the moment, the full NLO 
expression for including the 0{as) contributions to the matching conditions, 

is only available for the SM [13]. The contributions to Ek (i-e. to the imaginary part 
of the Hamiltonian) are also known for the two Higgs doublets model [14] and for 
the chargino contribution in the CMSSM [11]. In all these cases, the 0{as) terms in 
the matching amount to a very small correction, and the bulk of the effect is due to 
the running down to the hadronic scale. In our case, we find that the effect of NLO 
corrections to the running of the Left-Right operators is quite large (of the order of 
30% or more), which is somewhat worrying with respect to the convergence of the 
renormalisation-group improved perturbative expansion. 

The paper is organised as follows. Next section will introduce the in SUSY 

with the inclusion of the NLO QCD corrections in the evolution from the SUSY scale 
down to low energy. We will also provide the expression of the Wilson coefficients at the 
hadronic scale (to be detailed below) as a function of the Wilson coefficients and of as 
at the SUSY scale. Section 3 deals with the evaluation of the hadronic matrix elements 
of the local operators of Here we will replace the “traditional” values of the 

R-parameters in the VIA with the values which were recently obtained in a lattice 
computation [9]. Our quantitative results are presented in Sect. 4 and summarised in 
the Tables 2-7. These latter contain new constraints on the FC parameters S and a 
comparison with the previous results. The major effect of our improved computation 
is felt by Left-Right operators. As a further application, in Sect. 5, we consider the 
contribution to K^-K^ mixing in models with explicit R parity and lepton number 
violations. We obtain new bounds for the relevant R parity breaking A'-type couplings 
and compare them to those which were previously obtained in the VIA case and without 
QCD corrections. Finally Sect. 6 contains some conclusions and an outlook. 


2. Effective Hamiltonian for AS = 2 processes in SUSY 

In this Section, we describe the computation of gluino-mediated contributions to 
at the NLO in QCD. 

In ref. [7], was computed at the LO in three different cases: i) rrig ~ mg, 

ii) nig rrig and iii) nig 3> nig , where nig is the average squark mass and nig is the 
gluino mass. 

Case ii), in which rrig rrig, cannot be realized in the framework considered here, in 
which the soft SUSY breaking terms are introduced at the Planck scale. This is due to 
the fact that the evolution from the Planck to the electroweak scale forbids such a mass 
hierarchy. In fact, neglecting the effects of Yukawa couplings and A-terms, one obtains 
in the down sector the following approximate expression for the ratio x = at 


( 2 . 1 ) 


the electroweak scale, in terms of the value Xq at the Planck scale [15]: 

9x0 9 

^ ^ _ 

1 + 7xo 7 

Even if one starts at the superlarge scale with an extreme hierarchy between squark 
and gluino masses (xq 3> 1), at the electroweak scale the two masses will be of the 
same order. For this reason, we will not consider the case rriq rrig in our analysis. 
Case iii), in which rriq ^ rrig, can be realized in some special models, such as effective 
supersymmetry [16] or models with a light gluino [17]. However, due to the peculiar 
features and signatures of these models, a careful analysis of these cases would lie 
beyond the scope of this paper (see ref. [18] for a NLO analysis of AS = 2 processes 
in these models). Therefore, in the following we will consider the case rriq ~ rrig, 
and compare our results with the zeroth- [6] and leading-order [7] results previously 
published. 

We will perform our computation in the so-called mass insertion approximation [1]. 
One chooses the super-CKM basis for the fermion and sfermion states, where all the 
couplings of these particles to neutral gauginos are flavour diagonal, while the genuine 
SUSY FC effect is exhibited by the non-diagonality of the sfermion mass matrices. 
Denoting by the off-diagonal terms in the sfermion mass matrices (i.e. the mass 
terms relating sfermions of the same electric charge, but different flavour), the sfermion 
propagators can be expanded as a series in terms of 5 = A‘^/rfS, where rfi is the average 
sfermion mass. As long as is signihcantly smaller than ?Ti^, we can just take the 
hrst term of this expansion and, then, the experimental information concerning FCNC 
and CP violating phenomena translates into upper bounds on the Ss [4]-[6]. 

Obviously the above mass insertion method presents the major advantage that one 
does not need the full diagonalisation of the sfermion mass matrices to perform a test 
of the SUSY model under consideration in the FCNC sector. It is enough to compute 
ratios of the off-diagonal over the diagonal entries of the sfermion mass matrices and 
compare the results with the general bounds on the 6s which we provide here from 
AMk and Sk- This formulation of the mass insertion approximation in terms of the 
parameters 6 is particularly suitable for model-independent analyses, but involves two 
further assumptions: the smallness of the off-diagonal mass terms with respect to 
the diagonal ones, and the degeneracy of the diagonal mass terms in the super-CKM 
basis. The latter assumption, related to the use of the average squark mass rfi, is well 
justihed in our case, since, on quite general grounds, one does not expect a sizeable 
non-degeneracy of the first two generations of down-type squarks. It is also possible, 
however, to dehne the mass insertion approximation in a more general way, which is 
also valid for non-degenerate diagonal mass terms (see ref. [19] for details). 

There exist four different A mass-insertions connecting flavours d and s along a 
sfermion propagator: (^Af 2 ^ and (^Af 2 ^^^. The indices L and R 

refer to the helicity of the fermion partners. The amplitude for AS = 2 transitions in 
the full theory at the LO is given by the computation of the diagrams in Eg. 1. 




Figure 1: Feynman diagrams for AS = 2 transitions, with h,k,l,m = {L,R}. 

Having calculated the amplitude from the diagrams in £g. 1, one has to choose a 
basis of local operators and perform the matching of the full theory to the one described 
by We have adopted the form 

= E a g, + ^ a 4 (2.2) 

i=l i=\ 


where 


Qi = , 

Q2 = , 

Qs = dR^ldRSL , 

Q4 = , 

Q5 = , 


(2.3) 


and the operators gi, 2,3 are obtained from the gi, 2,3 by the exchange L ^ R. Here 
Qr,l = Pr,l d, with Pji L = (1 ± 75 ) 72 , and a and f3 are colour indices. 











At the lowest order in QCD, we obtain the following results for the Wilson coeffi¬ 
cients: 
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where x = and the functions fe{x) and fei^x) are given by: 


Mx) = 
hix) = 


6(1 + 3a;) Ina; + a;^ — 9a;^ — 9a; + 17 

6 (a; — 1 )^ 
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(2.4) 


(2.5) 


In the absence of the 0{as) corrections to the matching, we interpret the Ci given above 
as coefficients computed at the large energy scale Ms ~ nig ~ m^, i.e. Ci = Ci{Ms). 

The Next-to-Leading anomalous dimension matrix for the most general has 

been recently computed [8]. We use the Regularisation-Independent anomalous dimen¬ 
sion in the Landau gauge (LRI), since we will make use of matrix elements computed 
in lattice QCD with the same choice of renormalisation scheme (see ref. [8] for details 
on the computation). 

A full NLO computation would also require the 0{as) corrections to the matching 
conditions in eq. (2.4). Unfortunately, such corrections are not available yet. One might 
argue that, being of order as{Ms), these contributions should be small, as suggested 
by the cases of the SM and of the two Higgs doublet model; however, this statement 
can only be conhrmed by an explicit computation. Unfortunately, due to the absence 



of 0{as) corrections to the matching, our will be affected by a residual scheme 

dependence, which would be cancelled by the missing terms of order as{Ms). 

Eq. (2.4) is obtained by integrating out all SUSY 
particles simultaneously. We then have to evolve the 
coefficients down to the hadronic scale /r = 2 GeV, at 
which we have evaluated the matrix elements. The 
SM contribution can be computed independently and 
evolved from to fi using the well-known NLO 
QCD corrections [13]. 

We give here an analytic formula for the expres¬ 
sion of the Wilson coefficients at the scale /i = 2 GeV 
as a function of the initial conditions at the SUSY 

Table 1: Constants used in the ^his formula has been 

phenomenological analysis. obtained by using the values in Table 1 for the SM 

parameters. 

For Ms > rrit we obtain 

Crip) =T.Y. + iM) >(“• C,(Ms), ( 2 , 6 ) 

i s 

where, in the evolution of the coefficients from Ms, we have chosen Ms = {Mg + Mq)/2. 
rj = as{Ms)/Oisirrit), /i = 2 GeV and the magic numbers are given below:^ 

ai = (0.29, -0.69, 0.79, -1.1, 0.14) 

= (0.82, 0, 0, 0, 0), = (-0.016, 0, 0, 0, 0), 

6 p) = (0, 2.4, 0.011, 0, 0), = (0, -0.23, -0.002, 0, 0), 

^ _Q 03^ ^(23) ^ _o.018, 0.0049, 0, 0), 

^ ('Q^ -0.019, 0.028, 0, 0), = (0, 0.0028, -0.0093, 0, 0), (2.7) 

^ 0.0049, 0.43, 0, 0), = (0, 0.00021, 0.023, 0, 0), 

= (0, 0, 0,4.4, 0), = (0, 0, 0, -0.68, 0.0055), 

= (0, 0, 0,1.5, -0.17), = (0, 0, 0, -0.35, -0.0062), 

= (0, 0, 0, 0.18, 0), = (0, 0, 0, -0.026, -0.016), 

= (0, 0, 0, 0.061, 0.82), = (0, 0, 0, -0.013, 0.018), 

and we have only written the non-vanishing entries. The magic numbers for the evo¬ 
lution of Gi-a are the same as the ones for the evolution of Gi_ 3 . Formulae (2.6) and 
(2.7) can be used in connection with the 5-parameters evaluated at /i = 2 GeV, given 

^In the previously published version of this paper, JHEP 10 (1998) 008, the magic numbers in 
eq. (2.7) were erroneously given in the basis in eq. (13) of ref. [8] instead of the basis (2.3) used in 
this work (we thank P. Slavich and F. Zwirner for pointing this out to us). However, all the numerical 
results in Sect. 4 and 5 of the previous version were correct, since they were obtained using the full 
evolution matrix instead of eq. (2.6). 


Gonstants 


Values 


Ik 

md(2GeV) 

m,(2GeV) 


as{Mz 


497.67 MeV 
159.8 MeV 
7 MeV 
125 MeV 

1.3 GeV 

4.3 GeV 
175 GeV 

0.119 




in eq. (3.3), to compute the contribution to AM^ and Ek at the NLO in QCD for any 
model of new physics in which the new contributions with respect to the SM originate 
from extra heavy particles. One just has to plug in the expression for the Ci evaluated 
at the large energy scale Ms in his favourite model. When the 0{as) corrections to 
the Ci{Ms) are available, one can obtain a full NLO, regularisation-independent result; 
in the cases where this corrections have not been computed yet, the results contain a 
residual systematic uncertainty of order as{Ms). We note that, due to the presence 
of large entries in the NLO anomalous dimension matrix, a systematic uncertainty of 
a few percents is present in the QCD evolution from the SUSY scale to the hadronic 
one. For this reason we only present the signihcant hgures for the magic numbers. 


3. Hadronic Matrix Elements 


The matrix elements of the operators Qi between K mesons in the VIA are given by: 


{A-»|Qi|A'“>v,a = xV/f/I- , 


{A'»|Q2|A»)v,a = M ( 
{A»|Q3|A”)via = i ( 


M, 


K 


rus + rud 


MkPk, 


M, 


K 


TTls + rUd 


{K^\Q,\K\ia = 


(iL°|Q5|i^°)viA = 


M, 


K 


1 1 
24 4 KrUs + rUd 


21 


M, 


K 


1 1 
8 12 KrUg + rUd 


MkJk 1 


MkII, 


(3.1) 


where is the mass of the meson and rris and are the masses of s and d quarks 
respectively. Here and in the following, the same expressions of the H-parameters of 
the operators Qis are valid for the operators Q 1 - 3 , since strong interactions preserve 
parity. 

In the case of the renormalised operators, we dehne the H-parameters as follows: 


(^°l4(h)l^°) 


1 

24 \ms{ii) + md{^i) 


M, 
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24 
1 / 


ms{i^) +md(/i). 


M, 
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4 Vm^(/i) +md(/i)^ 
1 ( Mk 

12 I 


2 

MKfKB,{^^) 

2 

MKflB,{^i) , 


ms(/i) + md(/i) 


(3.2) 



where the notation (or simply Qi) denotes the operators renormalised at the 

scale /i. 

A few words of explanation are necessary at this point. The 5-parameter of the 
matrix element {K^\Qi\K^), commonly known as Bx, has been extensively studied on 
the lattice due to its phenomenological relevance [20], and used in many phenomeno¬ 
logical studies [21], For the other operators, instead, all the phenomenological analyses 
beyond the SM have used 5-parameters equal to one, which in some cases, as will be 
shown below, is a very crude approximation. 

In eq. (3.2) the operators and the quark masses are renormalised in the same scheme 
(RI, MS, etc.) at the scale fx and the numerical results for the 5-parameters, 
presented below refer to the Landau RI scheme. Moreover, without loss of generality, 
we have omitted terms which are of higher order in the chiral expansion, see eqs. (3.1), 
and which are usually included in the dehnition of the 5-parameters (see ref. [9] for a 
thorough discussion of the advantages of our dehnition of the 5-parameters). 

In our numerical study, we have used, for p = 2 GeV: 

5i(/i) = 0.60(6) (3.3) 

52(/i) = 0.66(4) 

BM = 1.05(12) 

54(/i) = 1.03(6) 

55(p) = 0.73(10) 


The central value we use for 5i = Bx corresponds to 5^‘^(2GeV) = 0.61, in 
agreement with the recent estimates of refs. [20, 22]. 52_5 have been taken from 

ref. [9], where all details of the computation can be found (for another determination 
of these 5-parameters, calculated with perturbative renormalization, see ref. [10]). 


4. Numerical Analysis of AS = 2 Processes 

We now present the results of a model-independent analysis of mixing, for the 

case rriq ~ rrig. The Kl~Ks mass difference AMx and the GP-violating parameter ex 
are given by: 


AMx = 2Re(A:°|5^^=2|R:°) , (4.1) 

The SUSY (gluino-mediated) contribution to the low-energy contains two 

real and four complex unknown parameters: rriq, rrig, ('^ 12 )('^ 12 )^^, ('^ 12 )^^ and 

{^ 12 ) These parameters can be determined once a specihc SUSY model is chosen, 
and the NLO contribution to AMx and ex can be directly computed by using the 



expression for the low-energy and the ^-parameters given in eqs. (2.4), (2.6) 

and (3.3). 

However, it is also useful to provide a set of model-independent constraints on the 
individual ^-parameters, obtained by neglecting the interference between the different 
SUSY contributions. This is justihed (a posteriori) by noting that the constraints on 
different ^-parameters in the Kaon case exhibit a hierarchical structure, and therefore 
interference effects between different contributions would require a large amount of hne 
tuning. 

The model-independent constraints are obtained by imposing that the sum of the 
SUSY contributions proportional to a single ^-parameter and of the SM contributions 
to AMk and Ek does not exceed the experimental value for these quantities. 

A comment on the SM contributions is necessary at this point. Let us hrst consider 
AMk- The short-distance SM contribution to this process is dominated by charm- 
quark exchange. Hence its value can be computed once the CKM elements Vcd and 
Vcs are determined. Making the very reasonable hypothesis that SUSY contributions 
to tree-level weak decays are negligible, the determination of the CKM elements Vcd 
and Vcs is unaffected by SUSY, and therefore the SM contribution to AMk can be 
computed even in the presence of large SUSY contributions to loop processes. We are 
not considering here long-distance contributions to AMk- The constraints we obtain 
from AMk are therefore more conservative if the long-distance SM contributions add 
up to the short-distance ones. 

On the other hand, the SM contribution to Ek depends on the phase in the CKM 
matrix. This phase is usually extracted, in the context of the SM, from the analysis of 
Ek and B^-B^ mixing. However, in the presence of large (unknown) SUSY contribu¬ 
tions to and B^-B^ mixing, the extraction of the CKM phase is not possible. 

We therefore treat it as a free parameter. Since the SM contribution to Ek is always 
positive, and vanishes for vanishing CKM phase, in order to obtain a conservative limit 
on SUSY contributions we set the CKM phase to zero, and allow the SUSY contribution 
to saturate the experimental value of Ek- 

We make use of the H-parameters given in eq. (3.3) and subtract one standard 
deviation to their central values, in order to extract a more conservative bound on 
SUSY parameters. 

Barring accidental cancellations, following the procedure explained above, we obtain 
the limits on the ^-parameters reported in Tables 2-7. We give results for different 
values of the average squark mass, since the naive scaling of the constraints like 1/m 
is modihed by the perturbative QCD corrections. 

For Left-Right mass insertions, we consider two possible (extreme) cases: | ( 5/2 )lr\ 
\i^i 2 )Rl\ and {Sf 2 )LR = {^ 12 )rl- In the second case, we combine the contributions pro¬ 
portional to (5i2)i_R, {^ 12 )RL and (5i2)Lij(<^i2)i?i- This approach improves the one of 
refs. [5] and [7], where contributions proportional to {Sf 2 ) lr{Si2)rl where considered 
independently from the ones proportional to (<5^2)iR and ((5f2)i?L- The interference 



between the various Left-Right terms, in the case ( 5 ^ 2 )^^? = produces a can¬ 

cellation effect which is particularly sizeable at the NLO for rriq = 500 GeV, around 
X = m?/m? = 1, as one can see from Tables 4 and 5. This means that our results are 
not reliable for this particular case, since interference effects can change dramatically 
the hnal results for small variations of the input parameters. 

On the other hand, the {Sf 2 )ll{Si 2 )RR contribution can be treated independently 
from the (5^2)iL ('^ 12 )??^ ones, since the latter ones do not generate Left-Right 

operators and are therefore suppressed with respect to the first contribution. 

Let us compare our results with the previous analyses of refs. [6] and [7]. In both 
papers, the matrix elements were computed in the VIA, whereas we are using the 
R-parameters computed in lattice QCD. This allows us a consistent matching of the 
renormalisation-scheme and p dependence of the QCD-corrected Moreover we 

keep into account important non-perturbative effects present in the hadronic matrix 
elements. 

If we use the same input parameters as in ref. [6], and neglect the SM contribution, 
we reproduce the results of ref. [6] for the QCD-uncorrected case. The values quoted in 
Tables 4-5 differ from the ones of ref. [6] due to the inclusion of the SM contributions 
in the present analysis and to the different choice of input parameters. As one can read 
from Table 4, our QCD-uncorrected and LO results are also different from those of 
ref. [7]. This is due to the different choice of the input parameters, to the neglect of the 
SM contribution in ref. [7] and to the choice of the scale fi. In particular, the results 
in Table 4 correspond to fi = 2 GeV, while the authors of ref. [7] chose to evolve the 
down to the scale fi' dehned by as{fi') = 1. This choice may be questionable, 
since perturbation theory is expected to break down at such low scales. 

5. Constraints on R-parity violating couplings 

In this section we update the constraints that can be derived from AS = 2 processes 
on some R-parity violating couplings. 

If R-parity invariance is not imposed, the superpotential may contain the following 
additional terms: 

Kjk''^id‘jdt. + X'ij^LiQjdl + XijkLiLjcl , (5.1) 

where i, j and k are generation indices, and we omitted possible bilinear R-parity 
violating terms [23]. The Xijk and are lepton number violating Yukawa couplings, 
while the X'Xj^ are baryon number violating ones. The Xijk are antisymmetric under the 
exchange of the hrst two generation indices, while the X'Xf. are antisymmetric under the 
interchange of the last two flavour indices. We consider the R-violating couplings in 
the super-GKM basis. 

Stringent constraints can be imposed on R-violating couplings by considering a 
variety of processes (see ref. [24] for a recent review on this subject). We concentrate 



on the limits that can be obtained from AS = 2 processes, and in particular on the 
ones arising from Left-Right effective operators. 

The most stringent constraints arise from the tree-level sneutrino-mediated contri¬ 
butions, that have the form 


Ci{mc) = 

i 




(5.2) 


Additional constraints can be obtained by considering box diagrams with four R- 
violating couplings. As was pointed out in ref. [25], however, stronger constraints on 
the product of two couplings can be obtained by considering box diagrams with the 
exchange of a slepton, together with a IR or a We consider only contributions 
with top-quark exchange, since when light quarks are present in the loop the OPE 
is modihed and the QCD evolution is much more involved. This computation is not 
available yet. In eq. (3) of ref. [25], the coefficient of the effective Hamiltonian at the 
electroweak scale was computed for mH± = tui, by integrating out simultaneously all 
heavy particles. In order to avoid the appearance of large logs of the form Inmf/m? in 
the matching conditions, we consider this contribution for mH± = rni = 200 GeV, and 
derive constraints on the product of the R-violating couplings 

Using the same procedure as that followed for R-conserving SUSY contributions, 
we obtain the constraints given in Table 8. As in the previous Section, we put the 
CKM phase to zero in order to obtain more conservative limits. 


6. Conclusions 

In our work we have provided an improved computation of the gluino-mediated SUSY 
contributions to mixing in the framework of the mass insertion method. The 

improvement consists in introducing the NLO QCD corrections to the [8] and in 

replacing the VIA R-parameters with their recent lattice computation [9]. As a glimpse 
at Tables 2-7 readily reveals, these improvements affect previous results in a different 
way, according to the operators of one considers. The effect is particularly large 

for Left-Right operators. 

We have also improved in the same way the constraints on some products of R- 
violating A'-type couplings that can be derived from mixing. Here the effects 

of QCD corrections, that were never considered before, are also very large, due to the 
Left-Right structure of the operator involved. 

We have provided an analytic formula for the most general low-energy at 

the NLO, in terms of the Wilson coefficients at the high energy scale. This formula can 
be readily used to compute AM^ and Ek in any extension of the SM with new heavy 
particles. 

The above results suggest several improvements for theoretical analyses. As for 77- 
mixing, a computation of the Wilson coefficients at Ms is needed to have the complete 



NLO QCD-corrected and scheme-independent result. Another issue that needs to be 
clarified is the role played by the other classes of FCNC contributions which are present 
in SUSY in addition to the gluino exchange. In particular the chargino-squark contri¬ 
butions should be considered, although, as we said, for generic squark mass matrices, 
we expect the gluino exchange to set the correct size of the whole FCNC contributions. 

Even more interesting is to extend this analysis to AB = 2 processes and to combine 
the limits which can be derived by considering FCNC processes involving B and K 
mesons simultaneously. The lattice derivation of the Bb parameters is in progress and 
so we are confident to be able to perform a similar analysis in the B case. 

FCNC and CP violating phenomena (in particular in B physics) are promising 
candidates for some indirect SUSY signal before LHC, and are in many ways comple¬ 
mentary to direct SUSY searches. From this point of view the theoretical effort to 
improve as much as possible our precision on FCNC computations in a SUSY model- 
independent framework is certainly worth and, hopefully, rewarding. 
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Table 2: Limits on Re A, B,C, D = (L, R), for an average squark 

rriq = 200 GeV and for different values of x = m‘~/m‘~. 
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X 
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Table 3: Limits on Im A, B,C,D = {L, R), for an average squark mass 

rriq = 200 GeV and for different values of x = 
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Table 4: Limits on Re i^ij)cD^ with A, B,C,D = {L, R), for an average squark mass 

rriq = 500 GeV and for different values of x = 
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Table 5: Limits on Im i^ij)cD^ with A, B,C,D = (L, R), for an average squark mass 

rriFi = 500 GeV and for different values of x = 
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Table 6: Limits on Re i^ij)cD^ with A, B,C, D = (L, R), for an average squark mass 

nipi = 1000 GeV and for different values of x = 
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Table 7: Limits on Im A, B,C,D = {L, R), for an average squark mass 

rriq = 1000 GeV and for different values of x = 



NO QGD, VIA 

LO, VIA 

LO, Lattice Bi 

NLO, Lattice Bi 

me, (GeV) 

ReA'2iA';2 

100 

1.4 X 10-30 

6.1 X 10-33 

6.3 X 10-33 

4.4 X 10-33 

200 

5.8 X 10-30 

2.2 X 10-30 

2.3 X 10-30 

1.6 X 10-30 

500 

3.6 X lO-o 

1.2 X lO-o 

1.3 X lO-o 

8.5 X 10-30 

mi (GeV) 

K 31 K 32 

200 

6.2 X 10-^ 

2.4 X 10-^ 

2.5 X 10-^ 

1.7 X 10-3 

me, (GeV) 


100 

3.1 X 10-32 

1.3 X 10-32 

1.4 X 10-32 

9.6 X 10-30 

200 

1.2 X 10-33 

4.8 X 10-32 

4.9 X 10-32 

3.4 X 10-32 

500 

7.8 X 10-33 

2.6 X 10-33 

2.7 X 10-33 

1.8 X 10-33 

mi (GeV) 

1 -1 

B 

to 

200 

1.3 X 10-3 

5.2 X lO-o 

5.4 X lO-o 

3.7 X 10-6 


Table 8 : Limits on \'i 2 iKi 2 on A' 3 ^A '32 from AMk and ek for different values of the 
relevant SUSY masses. 






































